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Abstract 

The failure time of samples of heterogeneous materials (wood, 
fiberglass) is studied as a function of the applied stress. It is shown 
that in these materials the failure time is predicted with a good ac- 
curacy by a model of microcrack nucleation proposed by Pomeau. It 
is also shown that the crack growth process presents critical features 
when the failure time is approached. 



PACS: 62.20.Mk , 46.30 



It is very well known that different materials subjected to a constant load 
may break after a certain time, which is a function of the applied load [1-5]. 
Many models have been proposed to predict this failure time, but the physical 
mechanisms remain often unclear |^. Very recently Pomeau proposed a 
model 1^, which explains quite well the failure time of microcrystals ^ and 
gels submitted to a constant stress. This model is based on the interesting 
idea that a nucleation process of microcracks has to take place inside the 
materials, in order to form the macroscopic crack. This nucleation process 
is controlled by an activation law, as the coalescence of a phase into another 
in a liquid-solid transition. Based on this prediction p|, L. Pauchard et al. 
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[0] found that if a constant load is applied to a bidimensional microcrystal, 
it breaks after a time r given by the equation r = Toe^°/^ , where P is 
the applied pressure, and Tq and Pq are constants. Bonn et al. |Q found a 
similar law for gels. Pomeau predicted that for three-dimensional microscopic 
systems the life-time should be 



r = rJ^^/Pr (1) 

where Tq is a characteristic time and Po a characteristic pressure, which 
mainly depend on the material characteristics, the experimental geometry 
and temperature. This idea is quite interesting and it merits to be checked 
experimentally in hetrogeneous materials, such as fiber glass and wood pan- 
nels. Indeed, in two recent papers [^, we have shown that in these 
materials the microcracks, preceding the main crack form something like a 
coalescence around the final path of the main crack. 

Driven by these observation we decided to study the behaviour of these 
materials as a function of time and to check whether eq.l could be useful 
in order to predict the sample life time. To do this we monitor the acoustic 
emission (AE) released before the final break-up by a sample placed between 
two chambers between which a pressure difference P is imposed. In Fig. 
1 a sketch of the apparatus is shown. We have prepared circular wood and 
fiberglass samples with a diameter of 22 cm. and a thickness of 4 mm (wood) 
and 2 mm (fiberglass). In our samples the AE consists of ultrasound bursts 
(events) produced by the formation of microcracks inside the sample. For 
each AE event, we record the energy e detected by the four microphones, 
the place where it was originated, the time at which the event was detected 
and the instantaneous pressure and deformation of the sample. The energy 
is defined as the integral of the sum of the squared signals. A more detailed 
description of the experimental methods can be found in P, p!0 . 



We first investigate the behaviour of the samples as a function of time 
when they are submitted to a constant load. Our interest is focused on the 
life-time of the sample, on the behaviour of the released acoustic energy near 
the fracture and on the distributions of energy and time elapsed between two 
consecutive events. The behaviour of the energy as a function of time for a 
system submitted to a constant load has been studied by geologists, but they 



were not specially interested in what happened near the fracture [11, ^ 



We first imposed a constant strain to our samples, as it has been made 
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for crystals. As strain is fixed, every microcrack leads to a pressure decrease, 
so the system reaches a stationary state. This is because a microfracture 
weakens the material. In the absence of microcraks the pressure remains 
constant. One sample was submitted to a large deformation (close to the 
fracture) and it did not break after three days. Therefore, at imposed strain, 
the effect observed in microcrystals is not valid for heterogeneous materials. 
On the other hand, if a constant stress is applied to the system, it will break 
after a certain time which depends on the value of the applied pressure. The 
reason for this is that after every single microcrack the same load must be 
endured by the weakened sample, so it becomes more and more unstable. 
We have submitted several samples to different constant pressures and we 
have measured the time until the break-up (the life-time r). The values 
obtained are well fitted by eq.(|l]), that is the exponential function predicted 
by Pomeau; the life-time expression r = ae~^^ proposed by Mogi |TT 



on 



the other hand, does not conform to our data[|T^. In Fig. 2a r is plotted 
versus in a semilog scale, and a straight line is obtained. Even if the 
pressure difference is very small the sample will eventually break, although 
the life-time can be extremely long. For example, using eq. (|l|) and the best 
fit parameters of fig. 2a, one estimate r ~ 5000 s at P = 0.43 atm. Halving 
the imposed pressure causes r to become extremely large : r = 4.4 ■ 10^^ 
years at P = 0.21atm). 

When a constant pressure is applied to the sample, the acoustic emission 
of the material is measured as a function of time. We find that the cumulated 
acoustic energy E diverg GS clS Qb function of the reduced time — , specifically 
E oc (^)'^ with 7 = 0.27 (see Fig. 2b). Notably, the exponent 7, found 
in this experiment with a constant applied pressure, is the same than the 
one corresponding to the case of constant stress rate Indeed it has been 
shown |10| that if a quasi-static constant pressure rate is imposed, that is 
P = Apt, the sample breaks at a critical pressure Pc and E realesed by the 
final crack precursors (microcrcaks) scales with the reduced pressure or time 
(time and pressure are proportional) in the following way: 

P — P T — f 
E (X {-^V = ( V (2) 

Pc T- 

where r = Pc/Ap in this case. Thus it seems that the real control param- 
eter of the failure process is time, regardless of the fact that either a constant 
pressure rate or a constant pressure is applied. 



3 



To find a general law, which is valid for a time dependent imposed stress, 
we intend to generalize the eq.|I| which is valid only for a constant imposed 
pressure. In the case where the pressure changes with time, it is reasonable 
to consider the entire history of the load. Therefore we consider that 

^exp(-(^r) 

' o 

is the density of demage per unit time. The certitude of breaking is obtained 
after a time r such that: 

r ie-(*)^rft = 1 (3) 

^0 To 

where and Pq have the previously determined value. Notice that this 
equation is equivalent to eq. ^ when a constant pressure is applied. 

To test this, we have applied the load to the sample following different 
schemes. We have first applied successive pressure plateaux in order to check 
whether memory effects exist. In fig. 3a the pressure applied to the sample 
is shown as a function of time. A constant load has been applied during a 
certain time ri, then the load is suppressed and then the same constant load 
for a time interval T2 is applied again. The sample breaks after a loading a 
time Ti + T2 which is equal to the time needed if the same load had been 
applied continuosly without the absence of load during a certain interval. 
Therefore there is a memory of the load history. The life-time formula (eq. 
^) is also valid if different constant loads are applied successively (fig. 3b). 
This concept can explain the violation of the Kaiser effect in these materials 

If the load is not constant, the life-times resulting from the proposed 
integral equation are still in good agreement with experimental data. A 
load linearly increasing at different rates Ap has been applied to different 
samples. The measured breaking times are plotted in fig.4 along with a 
curve showing the values computed from eq.^. Even if a quasi-static load is 
applied erratically (fig. 3d), the calculated life-time agrees with the measured 
one. These experiments show that eq. ^ describes well the life-time of the 
samples submitted to a time dependent pressure. 

The question is to understand why eq.(0) and (|]) works so well for a three 
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dimensional heterogeneous material. Indeed, in the Pomeau formulation 




Po = G\^\ (4) 



where Y is the Young modulus, T the temperature, K the Boltzmann constant 
and 7] the surface energy of the material under study. G is a geometrical 
factor which may depend on the experimental geometry, on defect shape and 
density. 

In our experiment, we found Pq = 0.62 atm for wood, which has Y=1.8 ■ 
10^ N/m^, and Po = 2.91 atm for fiberglass, which has Y=10^° N/m^. Thus 
the ratio between the values of Po found for the two materials is closed to 
the ratio of the square root of their Young modula. 

In contrast temperature does not seem to have a strong influence on r. 
In fact we changed temperature, from 300K to 380K which is a temperature 
range where the other parameters, Y and rj, do not change too much. For this 
temperature jump one would expect a change in r of of about 50% for the 
smallest pressure and of about 100% for the largest pressure. Looking at fig. 4 
we do not notice any change of r within experimental errors which are about 
10%. In order to maintain the change of r within 10% for a temperature 
jump of 80K one has to assume that the effective temperature of the system 
is about SOOOi^'. Notice that this claim is independent on the exact value of 
the other parameters and G. 

These observations seem to indicate that the nucleation process of mi- 
crocraks is activated by a noise much larger than the thermal one. Such a 
large noise can be probably produced by the internal random distribution of 
the defects in the heterogeneous materials that we used in our experiments. 
This internal random distribution of material defects evolves in time because 
of the appearance of new microcracks and the deformation of the sample. 
Therefore this internal and time dependent disorder of the material could 
actually be the mechanism that activates the microcrak coalescence and play 
the role of a very high temperature. Similar conclusions about a disordered 
induced high temperature have been reached in other disordered systems 



[ITSll . This is an important point that merits to be deeply explored. Simple 
numerical simulations which we performed in fuse networks seem to confirm 
this hypothesis. 

As a conclusion we have shown that a model based on the nucleation of 
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microcracks is in agreement with experimental data of failure of two het- 
erogeneous materials. This model seems to be quite general because it also 
explains the failure of gels |^ and microcrystals It will be certainly in- 
teresting to test it in other materials. However many questions remain open. 
The first one concerns the high temperature of the system. The second is 
related to the behaviour of acoustic energy close to the failure time. Indeed if 
AE is considered as a susceptibility |^ it is not easy to put together the ob- 
served critical divergency with a nucleation process. Probably the standard 
phase transition description can be only partially applied to failure because 
of the intrinsic irreversibility of the crack formation. All of these are very 
interesting aspects of this problem which certainly merit to be clarified in 
the future. 



References 

[1] Atkinson B. K., Introduction to fracture mechanics, Academic Press 
(1989). 

[2] Kaninnen F., Popelar C. H., Advanced Fracture Mechanics, Oxford Uni- 
versity Press, (1985). 

[3] Lawn B., Fracture of brittle Solids, Cambridge University Press(1993). 

[4] Shah S. P. Toughening Mechanisms in Quasi Brittle Materials, Kluwer 
Accademic Publisher (1991). 

[5] Dieter G., Mechanical Metallurgic, Mc Graw Hill, (1988). 

[6] Pomeau Y., Brisure spontanee de cristaux bidimensionnels courbes, C. 
R. Acad. Sci. Pans, vol. 314 II, 553-556 (1992), (in French). 

[7] Pauchard L. and Meunier J., Instantaneous and time-lag breaking of a 
two-dimensional solid rod under a bending stress, Phys. Rev. Lett. 70 
23, 3565-3568 (1993). 

[8] Bonn D., Kellay H., Prochnow M., Ben-Djemiaa K. and Meunier J., 
Delayed fracture of an inhomogeneous soft solid. Sciences, vol. 280, pp 
265-267 (1998). 



6 



[9] Garcimartin A., Guarino A., Bellon L. and Cilibcrto S., Statistical Prop- 
erties of fracture precursors, Phys. Rev. Lett. 79 17, 3202-3205 (1997). 

[10] Guarino A., Garcimartin A. and Cilibcrto S., An experimental test of 
the critical behaviour of fracture precursors, European Physical Journal 
B, 6, 13-24 (1998). 

[11] Mogi K., Study of elastic shocks caused by fracture of heterogeneous 
materials and its relation to earthquake phenomena. Bull. Earth Res. 
Inst. Tokyo 40, 125- (1962). 

[12] Stead D. and Szczepanik Z., Acoustic emission during uniaxial creep 
in potash. Fifth conference on acoustic emission/ micro seismic activity 
in geologic structures and materials, Pennsylvania State Univ., 97-111 
(1991). 

[13] Maes C. Van Moffaert A., Frederix H. and Strauven H., Criticality 
in creep experiments on cellurar glass, Phys. Rev. B 57 9, 4987-4990 
(1998). 

[14] Several functions with many free parameters have been proposed [1-5] 
in order to fit the dependence of the sample life time on the applied 
stress. Among the fitting functions, which are supported by a micro- 
scopic physical mechanism we show in fig. 2a the two functions which 
give the best fit of our data. 

[15] Sollich P., Rhcological constitutive equation for a model of soft glassy 
materials, Phys. Rev. E, 58, 1 738-759 (1998). 

[16] Zapperi S., Vespignani A. and Stanley H.E., Plasticity and avalanche be- 
haviour in microfracturing phenomena. Nature, vol 338, 658-660 (1997). 

Correspondence should be addressed to S. Cilibcrto (e-mail: cilibe@physique.ens- 
lyon.fr) 



7 



Figure Captions 



1. Sketch of the apparatus. S is the sample, DS is the inductive displace- 
ment sensor (which has a sensitivity of the order of 1 fim). M are 
the four wide-band piezoelectric microphones. P=Pi-P2 is the pressure 
supported by the sample. P is measured by a differential pressure sen- 
sor ( sensitivity = 0.002 atm) that is not represented here. EV is the 
electronic valve which controls P via the feedback control system Ctrl. 
HPR is the high-pressure air reservoir. 

2. Measures on wood samples: a) The time r needed to break the wood 
samples under an imposed constant pressure P is here plotted as a 
function of 1/P'^ in a semilog scale. The dashed line represents the 
solution proposed by Mogi [|n|](r = ae~^^). The continuous line is the 
solution proposed by Pomeau for microcrystals (r = TqC^^"^^^ ). In 
the plot To = 50.5 s and Po = 0.63 atm. Every point is the average 
of 10 samples. The error bar is the statistical uncertainty. For the 
fiberglass samples, we find To = 44.6 s and Po = 2.91 atm. b) The 
cumulated energy E, normalized to Emax, as a function of the reduced 
control parameter ^ at the neighborhood of the fracture point (Case 
of imposed constant pressure). The circles are the average for 9 wood 

samples. The solid line is the fit E = Eq ^- The exponent found, 

7 = 0.26, does not depend on the value of the imposed pressure. In the 
case of a constant pressure rate the same law has been found [0, ^. 

3. The imposed time dependent pressure (bold dotted line) is plotted as 
a function of time in the case of wood samples. The continuous line is 
the integral in time of the function / (P) = ^e~^°/^ . On the basis of 
eq. ^ the predicted breaking time r is obtained when the integral of 
/ (P) is equal to 1. The horizontal distance between the two vertical 
dashed lines in each plot represent the difference between the predicted 
and the measured breaking time. In a) a constant pressure has been 
applied during about 700 s, then the load is suppressed and then the 
same constant load is apphed again. The difference between the life- 
time predicted by (eq. |^) and the experimental result is of 3%. b) 
Here two pressure plateaux of different value are successively applied 
to the sample. The difference between the measured and the predicted 
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life-time is of 5%. In c) an erratic pressure is applied to the sample. 
Here the error is of 10%. 

4. A load linearly increasing at different rates Ap has been applied to dif- 
ferent samples. The measured breaking times are plotted as a function 
of Ap in a loglog scale; circles and squares represent the measures on 
wood and fiberglass samples respectively at T=300 K. Bold triangles 
represent measures on wood samples at T=380 K. The lines are the 
life time calculated from eq.|^ using the best fit values for Po and Tq. 
These experiments show that eq. ^ describes well the life-time of the 
samples submitted to a time dependent pressure. 
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